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USERSIN 


AN  INTERACTIVE  USER- I NTERFACE  FOR 


FORTRAN  SIN 


1 .  INTRODUCTION 


An  increasing  requirement  for  hydrodynamic  calculations  to  predict 
explosive  behaviour  and  weapon  performance  has  arisen  from  the  complexity  and 
expense  of  detailed  experimentation  in  these  fields.  Many  explosive  and 
weapon  systems  can  be  simulated  by  relatively  simple  one-dimensional  computer 
codes.  Phenomena  such  as  explosive -metal  interactions,  fragmentation  and 
shock  studies  are  ideally  suited  to  this  approach. 


For  some  years,  MRL  has  had  access  to  SIN  ,  a  Lagrangian  one¬ 
dimensional  hydrodynamic  program,  used  extensively  overseas  for  problems- 
including  chemical  reactions,  material  failure  and  phase  transitions.  Fortran 
SIN  was  developed  at  the  Los  Alamos  Scientific  Laboratories  in  1964  for  a 
CDC6600  computer  and  has  been  modified  to  run  on  the  CSIRO  computer  network 
which  is  based  on  a  Cyber  7600. 


Running  the  SIN  code  directly  on  the  Cyber  76  clearly 
demonstrated  the  inherent  difficulties  in  compiling  the  input  data  for  a  batch 
run.  This  Report  then  describes  USERSIN  (Fig.  1),  a  program  developed  to 
create  interactively  the  input  data  deck  for  SIN. 


2. _ SIN: 


The  SIN  code  requires  a  large  amount  of  input  data  defining  material 
properties,  the  number  of  cells  over  which  the  calculations  are  to  be  made,  as 
well  as  the  time  increment  to  be  considered.  Further  input  data  include  the 
parameters  for  elastic-plastic  calculations  on  the  solids  and  those  parameters 
defining  the  equation  of  state  of  the  detonation  products,  the  latter 
constants  being  evaluated  by  the  BKW  code'3'.  (Appendix  1  lists  the  input 
data  required  for  SIN). 


As  it  is  necessary  to  collect  and  format  correctly  all  the  data 
before  SIN  can  be  run,  such  a  procedure  is  both  a  time-consuming  and  error- 
prone  operation. 


3.  USERSIN: 


The  user  interface  usersin  written  in  extended  Fortran  4  IV. 
Usersin  takes  advantage  of  the  computer's  connected  input/output  subsystem 
(CIO)  to  create  interactively  the  data  deck  for  subsequent  use  by  SIN. 


The  program  is  menu  driven.  This  means  that  it  displays  a  table  of 
options,  to  which  a  user  is  able  to  select  any  option  by  entering  the 
corresponding  key.  Hie  user  moves  from  one  menu  to  the  next  being  promoted 
with  a  suitable  message  to  enter  the  required  data.  These  data  are  then 
validated  and  written  into  the  data  deck.  The  program  also  assumes  default 
values  which  may  be  changed  by  the  user. 


The  equation  of  state  parameters  and  other  physical  material 
constants  are  stored  in  a  data  base  file.  This  file  is  a  sequential  file  with 
the  name  of  the  material  its  key  (i.e.  TNT).  The  file  is  automatically 
searched  after  the  components  name  is  entered.  If  the  component  is  on  file, 
the  appropriate  parameters  are  retrived  and  placed  into  the  data  deck, 
otherwise  the  user  is  solicited  for  these  parameters. 


Users  of  this  program  are  requested  to  store  the  components  they  use 
into  this  data  base  file  when  prompted,  thus  increasing  the  usefulness  of  this 
data  base  to  future  users. 


The  Usersin  package  is  simple  to  master,  making  SIN  a  more  reliable 
and  efficient  research  tool. 


A  typical  run  of  Usersin  for  an  expanding  steel  casing  is  shown  in 
Appendix  3. 


4 .  MAINTENANCE 


For  ease  of  maintaining  USERSIN  the  program  has  been  fully 
commented,  each  subroutine  having  one  particular  function  which  is  explained 
in  detail  at  the  routines  beginning.  The  program  was  written  to  be  as  machine 
independent  as  possible  and  a  comment  appears  where  it  was  not  possible.  The 
input/output  variable  names  have  been  kept  the  same  as  those  used  in  SIN  so 
the  maintenance  programmer  can  easily  swap  from  one  to  the  other. 
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Hie  data  base  file  is  under  the  format  of 
( 1  HZ,  2A6,  2X  15,/,  2 (1 X,  E18.11),  /,  9(4(1X,  E18.11)/) 


KEY:-  NAME  2A6 

EXPLOSIVE  OR  SOLID  (1  OR  0)  15 

DENSITY  El  8.11 

MASS  FRACTION  El  8.11 

C,S  Cl,  SI,  CV  SPALLP,  USP,  SPA,  F,  G,  H,  I,  J,  Mo,  2/3  YO, 

Y,  MINV,  VSVJ,  PLAP,  A,  B,  C,  D,  E,  K,  L,  M,  N,  CV1  ,  Q,  R,  S,  T, 

U,  Z0,  VFACT,  EO,  Z  El  8.11 


5.  CONCLUSION 


Development  of  USERSIN  has  been  completely  justified  by  its  ease  of 
operation,  its  data-storage  capability  and  its  use. 
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USE  OF  SIN 


APPEND  I  X  1 


DATA  REQUIRED  FOR  SIN 


Variable  Default 

name  value 


Units 


Def In  It  Ions 


TITl  Time  -  date 


A  suitable  title  for  the  run  (0-80  char) 


IHC 

No 

** 

ALPHA 

SUB 

IPINC  400 

IGINC  200 

I  LB  00N 


Is  a  heat  calculation  required  IY/N) 

The  geometry  of  the  system 
SLAB  -  slab  geometry 
CfL  -  cylindrical  geometry 
SPHER  -  spherical  geometry 

Every  IPINC  cycles  are  printed 

Every  IGINC  cycles  are  graphed 

Left  hand  boundary  conditions 

CON  -  Contlnuun 

FRE  -  free  space 

INIT  -  Initial-final  piston 

A+BT  -  A  +  BT  piston 

STDY  -  steady  state  reaction  zone 
piston 


IRB  OON  -  Right  hand  boundary  conditions 

CON  -  contlnuun 

FRE  -  free  space 

INIT  -  Initial-final  piston 

A+BT  -  A  +  BT  piston 

Only  required  If  ILB,  IRB  Is  equal  to  INIT  "Initial-final  piston" 

ALB,LRB  -  cmAis  Initial  piston  velocity 


BIS,  BRB 


cm/ps 


final  piston  velocity 


Variable 

name 


Default 

value 


Units 


Definitions 


Only  required  If  I  LB,  IRB  Is  equal  to  AtBT  "A  +  BT  piston" 


ALB,  ARB 

— 

• 

A 

BLB,  BRB 

- 

- 

B 

On  1  y  req  u  1  red 

If  1  LB ,  IRB  1  s  equal 

to  1  NIT,  A+BT  or  STDY 

NLI,  NRI 

- 

- 

The  nunber  of  cel  Is  Initial  ly  In  the 
calculation  from  the  left,  right 

If  0  al  1  cel  Is  are  taken. 

NLINC,  NR  INC 

- 

- 

Every  NLINC,  NRINC  cycles  a  cel  1 

Is  added  on  from  the  left,  right 

The  fol lowing 

data  Is  required  tor 

al 1  components 

NAM 

-  X 

- 

Name  of  the  component  (0-12  char) 

NINC 

-  X 

- 

Nunber  of  space  Increments  for  this 
component 

DX 

- 

cm 

Size  of  these  cells 

DTIME 

0.2  of  DX 

US 

Size  of  the  time  Increnents 

IEXP 

-  X 

- 

Is  the  component  a  gas  or  explosive 

1 Y/NJ? 

1  VIS 

PIC 

PIC  -  Pic  viscosity 

LAND  -  Lands  hot  f  form 

REAL  -  Real  form 

VFACT 

2 

- 

Viscosity  factor 

PO 

1  x  10-6 

MBAR 

Initial  presstre 

TO 

300 

K 

Initial  temperature 

EO 

0 

MBar-cq/g 

Initial  Internal  energy 

Variable 

name 

Default 

value 

Units 

Definitions 

uo 

0 

cm/ys 

Initial  velocity 

RHOO 

(KRHOOOOX 

g /cm3 

Initial  density 

WO 

(XWIX1  X 

Initial  mass  traction 

VSW 

0.01 

Solid  paraneters 

C,S 

-  X 

- 

Co-ef f Iclents  to  a  linear  set  of  Us 
and  Up  used  from  the  Initial  pressure 
to  the  switch  pressure. 

Cl,  SI 

-  X 

' 

Second  set  of  co-ef f  Iclents  to  the 
linear  fit  of  Us  and  Up  used  from  the 
switch  pressure  to  the  max.  pressure 

cv+ 

X 

Ca  l/g/ Degree 

Heat  of  capacity  of  the  component 

VO 

l/DENS 

cm3/g 

Initial  volume  of  the  condensed 
component 

ALPH 

-  X 

- 

Linear  co-ef flclent  of  thermal 
expansion 

F.G.H.J.I, 

-  X 

- 

Qonstants  used  1 n  the  Horn  condensed 
components  equation 

GAMMA 

-  X 

- 

■y-law  gas  constant 

Not  requested 

If  the  component  Is  a 

gas  or  an  explosive 

SPALLP 

-  X 

MBar 

Interface  spalling  pressure 

USP 

-  X 

MBar 

Ultimate  spalling  pressure 

SPA 

-  X 

- 

Spal ling  constant  to  relate  pressure 
and  tension  rate 

SHEAR 

-  X 

MBBr 

Shear  modulus 

TO 

-  X 

MBar 

Yield  strength 

Only  requested  If  tit*  component  Is  a  gas  or  explosive 


* 


Variable 

name 

Default 

value 

Units 

Definitions 

z 

1  El  5  X 

Frequency  factor  for  an  Arrhenius  birn 

E 

4.5E4  X 

cal/mole 

Activation  energy  for  an  Arrhenius  burn 

VCJ 

X 

cm^/g  or  on/ps 

Volume  of  the  component  for  a  C-J  burn 
or  the  C-J  detonation  velocity  for  a 
Ganma-Law  Taylor  wave. 

A#8fC  *  D ,  E 

-  X 

Parameters  defining  the  equation  of 

K,L,M,N,0 

-  X 

- 

state  of  the  detonation  products. 

0,R.S,T,U 

-  X 

Parameters  obtained  from  the  BKW  code^ 

Z1 

-  X 

- 

A  constant  used  to  change  the  standard 
state  to  be  consistent  with  the  solid 
explosive  standard  state  usual  ly  s  0.1 

IBRN 

ARR  X 

ARR  -  Arrhenius  burn 

C-J  -  C-J  volume  burn 

GAMMA  -  Gamma-Law  Taylor  wave  burn 

FOREST  -  Forest  Fire  burn  (when  available) 

Al  1  variables 

wl+h  an  "X"  marked 

In  the  default  value 

column  are  store<  on  the  data 

base  file. 


issed  elsewhere  in  this  report  a  user  interface  has  been 
the  data  deck  for  the  SIN  code. 


Jed  binary  version  of  USERSIN  are  stored  on  the  library 
inly  necessary  to  execute  this  file  to  run  USERSIN.  This  must 
an  interactive  subsystem  of  the  Cyber,  once  logged  on. 


ET  (USERSIN) 


|tnd  a  brief  description  of  the  program  should  then  be 


Questions  asked  and  at  the  end  of  the  session  the  data 
cards  will  have  been  created  and  will  only  require 


jin  is  a  menu  driven  program,  for  example  the 
fcng  initial  menu.  The  user  simply  enters  the 


quit  edit  subsystem 
gain  entry  to  CIO 
excute  USERSIN 


quit  CIO  subsystem 

gain  entry  to  EDIT 

copy  the  data  deck  into  work  space 

excute  SIN 


DEFAULT 


ter  title 
tions 

the  system 

dary  conditions 
dary  conditions 


Date  and  Time 
Not  included 
SLAB 

CONTINUUM 

CONTINUUM 


option  at  least  the  first  three  characters 
rent  value  of  the  option  is  given  on  the 


tered  the  following  options  are  offered. 


OWING: 


ETRY 
ICAL  GEOMETRY 
GEOMETRY 


Ihe  program  will 

validate  the  reply  and  return  to  the 

initial  menu 

changing  the  default  value 

as 

requested. 

INITIAL  MENU 

DEFAULT 

TITL 

- 

80  char  title 

Date  and  Time 

HEAT 

- 

Heat  calculation 

Not  included 

GEOM 

- 

Geometry  of  the  system 

CYLINDRICAL 

LHB 

- 

Left  hand  boundary  conditions 

CONTINUUM 

RHB 

- 

Right  hand  boundary  conditions 

CONTINUUM 

CONT 

- 

to  continue 

To  move  to  the  next  menu  CON  (continue)  must  be  entered.  The  user 
works  his  way  through  the  menu's  specifying  each  parameter,  Usersin  reads  and 
validates  his  reply  before  creating  the  data  deck. 


APPENDIX  3 


System  ? 

:  cio 

Cio  version  15/10/82 
: submit, usersin 
Enterins  7000  mode... wait 

THIS  IS  THE  FRONT  END  PACKAGE  FOR  SIN. 

A  ONE-DIMENSIONAL  HYDRODYNAMIC  CODE  FOR  PROBLEMS 
WHICH  INCLUDE  CHEMICAL  REACTIONS,  ELASTIC-PLASTIC  FLOW 
, SPALLING,  AND  PHASE  TRANSITIONS. 

IT  IS  A  QUESTION  ANSWER  TYPE  PROGRAM,  CONTROL  'H'  WORKS  AS  BACKSPACE 
IF  THE  WRONG  DATA  IS  ENTERED  EITHER  AN  ERROR  MESSAGE  WILL  APPEAR  OR  YOU  CAN 
EDIT  THE  ENTERED  DATA  AT  THE  END  OF  EACH  SECTION 


DEFAULT 

CURRENTLY 

NOT  INCLUDED 
SLAB  GEOMETRY 
A  CONTINUUM 
A  CONTINUUM 


s  tit 


ENTER  A  TITLE  FOR  THIS  RUN  s 

.-cylinder  test  comp  b  steel  4.225,  5.225  cm  dimensions 


ENTER  ONE  OF  THE  FOLLOWING 

TITLE  -  FOR  THE  TITLE  OF  THIS  RUN  (0-80  CHAR) 
04/03/83  16.58.48. 

HEAT  -  IF  HEAT  CALCULATIONS  ARE  REQUIRED 
GEOM  -  FOR  THE  GEOMETRY 
LHB  -  LEFT  HAND  BOUNDARY  CONDITIONS 
RHB  -  RIGHT  HAND  BOUNDARY  CONDITIONS 
CONT  -  TO  CONTINUE 


ENTER  ONE  OF  THE  FOLLOWING 


DEFAULT 


TITLE  -  FOR  THE  TITLE  OF  THIS  RUN  (0-80  CHAR) 
CYLINDER  TEST  COMP  B  STEEL  4.225,  5.225  CM  DIMENSIONS 
HEAT  -  IF  HEAT  CALCULATIONS  ARE  REQUIRED 
GEOM  -  FOR  THE  GEOMETRY 
LHB  -  LEFT  HAND  BOUNDARY  CONDITIONS 
RHB  -  RIGHT  HAND  BOUNDARY  CONDITIONS 
CONT  -  TO  CONTINUE 


CURRENTLY 

NOT  INCLUDED 
SLAB  GEOMETRY 
A  CONTINUUM 
A  CONTINUUM 


:geo 


ENTER  GEOMETRY 


SLAB  -  FOR  SLAB  GEOMETRY 
CYLIND  -  FOR  CYLINDRICAL  GEOMETRY 
SPHER  -  FOR  SPHERICAL  GEOMETRY 

:cyl 


ENTER  ONE  OF  THE  FOLLOWING 


DEFAULT 


TITLE  -  FOR  THE  TITLE  OF  THIS  RUN  (0-80  CHAR) 
CYLINDER  TEST  COMP  B  STEEL  4.225,  5.225  CM  DIMENSIONS 
HEAT  -  IF  HEAT  CALCULATIONS  ARE  REQUIRED 
GEOM  -  FOR  THE  GEOMETRY 
LHB  -  LEFT  HAND  BOUNDARY  CONDITIONS 
RHB  -  RIGHT  HAND  BOUNDARY  CONDITIONS 
CONT  -  TO  CONTINUE 


CURRENTLY 

NOT  INCLUDED 
CYLINDRICAL  GEOMETRY 
A  CONTINUUM 
A  CONTINUUM 


:cont 


OUTPUT  MENU 


ENTER  ONE  OF  THE  FOLLOWING 


DEFAULT 


: print 


PRINT 
GRAPH 
DUMP 
IPR  - 
ITR  - 
IWR  - 
HR  - 
I  UR  - 

CONT 


-  CYCLES  PRINTED  400 

-  CYCLES  GRAPHED  .400E+03 

-  CYCLES  DUMPED  .200E+07 

FOR  A  PLOT  OF  PRESS  VS  RADIUS 

FOR  A  PLOT  OF  TEMP  VS  RADIUS 
FOR  A  PLOT  OF  MASS  FRAC  VS  RADIUS 
FOR  A  PLOT  OF  VOLUME  VS  RADIUS 
FOR  A  PLOT  OF  PART  VEL  VS  RADIUS 

NO  GRAPHICAL  OUTPUT  REQUESTED 


-  TO  CONTINUE 


ENTER  N,  WHERE  EVERY  N  CYCLES  ARE  PRINTED 
slOO 


OUTPUT  MENU 


ENTER  ONE  OF  THE  FOLLOWING 


DEFAULT 


PRINT  -  CYCLES  PRINTED 

GRAPH  -  CYCLES  GRAPHED 

DUMP  -  CYCLES  DUMPED 

IPR  -  FOR  A  PLOT  OF  PRESS  VS  RADIUS 

ITR  -  FOR  A  PLOT  OF  TEMP  VS  RADIUS 

IWR  -  FOR  A  PLOT  OF  MASS  FRAC  VS  RADIUS 

HR  -  FOR  A  PLOT  OF  VOLUME  VS  RADIUS 

I UR  -  FOR  A  PLOT  OF  PART  VEL  VS  RADIUS 


100 

.400E+03 

.200E+07 


NO  GRAPHICAL  OUTPUT  REQUESTED 


-  TO  CONTINUE 


OUTPUT  MENU 


ENTER  ONE  OF  THE  FOLLOWING 


DEFAULT 


PRINT  -  CYCLES  PRINTED 


100 

GRAPH  -  CYCLES  GRAPHED  .400E+03 

DUMP  -  CYCLES  DUMPED  .200E+07 

IPR  -  FOR  A  PLOT  OF  PRESS  VS  RADIUS 

ITR  -  FOR  A  PLOT  OF  TEMP  VS  RADIUS 

IWR  -  FOR  A  PLOT  OF  MASS  FRAC  VS  RADIUS 

IIR  -  FOR  A  PLOT  OF  VOLUME  VS  RADIUS 

I UR  -  FOR  A  PLOT  OF  PART  VEL  VS  RADIUS 

A  GRAPH  OF  PART  VEL  VS  RADIUS 
WILL  BE  PLOTTED  EVERY  .400E+03  CYCLES 

CONT  -  TO  CONTINUE 


ENTER  THE  NAME  OF  COMPONENT  1 ,  OR 
OF  THE  COMPONENTS  ON  FILE 
s? 

THE  COMPONENTS  ON  FILE  ARE 

COMPONENTS  ON  FILE 

STEEL 

RDX 

TNT 

H6 

COMP  B 

AIR 

PETN 

POLYETHYLENE 

WATER 

PETN(RHO-1 .4 
PENTOLITE 


A  QUESTION  MARK  ’? '  FOR  A  LIST 


COMMENT 

PETERS  PROBLEM  RDX  FILLED 
PETERS  PROBLEM  RDX  FILLED 
INITIATION  OF  TNT  (WATER  ONLY) 
MARK  82  H6  FULL  SIZE 
4.5  INCH  SHELL/BARREL/ JACKET 

DAVID’S  PROB 

DAV'S 

DAV'S 

DAV’S 

DAV’S 


I 


THERE  ARE  11  COMPONENTS  ON  THE  COMPONENTS  PILE 

ENTER  THE  NAME  OF  COMPONENT  1 
:  corap  b 

COMP  B  HAS  BEEN  FOUND  ON  THE  COMPONENTS  FILE 

ENTER  THE  NUMBER  OF  SPACE  INCREMENTS  FOR  COMP  B 
:50 


ENTER  THE  SPACE  INCREMENT  SIZE  FOR  COMP  B  (CM] 

-.0.0845 


A  GOOD  APPROXIMATION  OF  THE  TIME  INCREMENT  IS  0.2  OF  THE  CELL  WIDTH 
WHICH  IN  THIS  CASE  IS  .1690000E-01  IS  THIS  APPROXIMATION 

VALID  [Y/N] 

sy 


PLEASE  ENTER  ONE  OF  THE  FOLLOWING 


DEFAULT 


vise 

- 

FOR  THE  VISCOSITY  FACTOR 

2.00 

BURN 

- 

FOR  TYPE  OF  BURN 

ARRHENIUS  1 

BURN 

FORM 

- 

FOR  TYPE  OF  VISCOSITY 

PIC  VISCOSITY 

PRES 

- 

INITIAL  PRESSURE 

.1  00E-05 

MBAR 

TEMP 

- 

INITIAL  TEMPERATURE 

.300B4-03 

K 

ENGY 

- 

INITIAL  INTERNAL  ENERGY 

0. 

MBAR-CC/G 

VELC 

- 

INITIAL  VELOCITY 

0. 

CM/MICROSEC 

FREQ 

- 

FREQUENCY  FACTOR 

.100E+16 

ACTV 

- 

ACTIVATION  ENERGY 

.450E+05 

NUMB 

- 

NUMBER  OF  CELLS 

50 

SIZE 

- 

SIZE  OF  THE  CELLS 

.845B-01 

CM 

TIME 

- 

TIME  INCREMENTS 

.169E-01 

MICROSEC 

MASS 

- 

MASS  FRACTION 

0. 

G 

DENS 

- 

DENSITY 

.1  65E+01 

G/CC 

PAR 

CONT 

: 

PARAMETERS  RELATING  TO  THE 
TO  CONTINUE 

BQU  OF  STATE 

:cont 


ARE  THERE  ANY  MORE  COMPONENTS  ?  (Y/N] 

:y 


ENTER  THE  NAME  OF  COMPONENT  2,  OR  A  QUESTION  MARK  FOR  A  LIST 
OF  THE  COMPONENTS  ON  FILE 

t steel 


STEEL 


HAS  BEEN  FOUND  ON  THE  COMPONENTS  FILE 


ENTER  THE  NUMBER  OF  SPACE  INCREMENTS  FOR  STEEL 
:40 


ENTER  THE  SPACE  INCREMENT  SIZE  FOR  STEEL  [CM] 

:0.025 

A  GOOD  APPROXIMATION  OF  THE  TIME  INCREMENT  IS  0.2  OF  THE  CELL  WIDTH 
WHICH  IN  THIS  CASE  IS  . 500000 OE-02  IS  THIS  APPROXIMATION 

VALID  [Y/N] 

:y 


PLEASE  ENTER  ONE  OF  THE  FOLLOWING 

DEFAULT 


vise 

- 

FOR  THE  VISCOSITY  FACTOR 

2.00 

FORM 

- 

FOR  TYPE  OF  VISCOSITY 

PIC  VISCOSITY 

PRES 

- 

INITIAL  PRESSURE 

.1 00E-05 

MBAR 

TEMP 

- 

INITIAL  TEMPERATURE 

.300E+03 

K 

ENGY 

- 

INITIAL  INTERNAL  ENERGY 

0. 

MBAR-CC/G 

VELC 

- 

INITIAL  VELOCITY 

0. 

CM/MICROSEC 

NUMB 

- 

NUMBER  OF  CELLS 

40 

SIZE 

- 

SIZE  OF  THE  CELLS 

•250E-01 

CM 

TIME 

- 

TIME  INCREMENTS 

.500E-02 

MICROSEC 

MASS 

- 

MASS  FRACTION 

.1 00E+01 

G 

DENS 

- 

DENSITY 

.792E4-01 

G/CC 

PAR 

- 

PARAMETERS  RELATING  TO  THE  EQU  OF 

STATE 

CONT 

- 

TO  CONTINUE 

:cont 


ARE  THERE  ANY  MORE  COMPONENTS  ?  [Y/N] 

:y 


ENTER  THE  NAME  OF  COMPONENT  3,  OR  A  QUESTION  MARK  '?’  FOR  A  LIST 

OF  THE  COMPONENTS  ON  FILE 

:air 


AIR  HAS  BEEN  FOUND  ON  THE  COMPONENTS  FILE 


ENTER  THE  NUMBER  OF  SPACE  INCREMENTS  FOR  AIR 
:10 


ENTER  THE  SPACE  INCREMENT  SIZE  FOR  AIR  [CM] 

:1  .0 

A  GOOD  APPROXIMATION  OF  THE  TIME  INCREMENT  IS  0.2  OF  THE  CELL  WIDTH 
WHICH  IN  THIS  CASE  IS  .2000000E+00  IS  THIS  APPROXIMATION 

VALID  [Y/N] 

:y 


7 


PLEASE  ENTER  ONE  OF  THE  FOLLOWING 

DEFAULT 


vise 

- 

FOR  THE  VISCOSITY  FACTOR 

2.00 

BURN 

- 

FOR  TYPE  OF  BURN 

ARRHENIUS  BURN 

FORM 

- 

FOR  TYPE  OF  VISCOSITY 

PIC  VISCOSITY 

PRES 

- 

INITIAL  PRESSURE 

.1  00E-05 

MBAR 

TEMP 

- 

INITIAL  TEMPERATURE 

.300B+03 

K 

ENGY 

- 

INITIAL  INTERNAL  ENERGY 

0. 

MBAR-CC/G 

VELC 

- 

INITIAL  VELOCITY 

0. 

CM/MICROSEC 

FREQ 

- 

FREQUENCY  FACTOR 

0. 

ACTV 

- 

ACTIVATION  ENERGY 

0. 

NUMB 

- 

NUMBER  OF  CELLS 

10 

SIZE 

- 

SIZE  OF  THE  CELLS 

.100E+01 

CM 

TIME 

- 

TIME  INCREMENTS 

.200E+00 

MICROSEC 

MASS 

- 

MASS  FRACTION 

.1 00E4-01 

G 

DENS 

- 

DENSITY 

.1 29E-02 

G/CC 

PAR 

CONT 

— 

PARAMETERS  RELATING  TO  THE 
TO  CONTINUE 

EQU  OF  STATE 

:COnt 


ARE  THERE  ANY  MORE  COMPONENTS  ?  [Y/N] 
sn 

DO  YOU  WISH  TO  PLACE  ANY  OF  THE  COMPONENTS  ON  THE  FILE  [Y/N] 

##  REPLY  SHOULD  BE  'Y'  OR  'N'  ##  TRY  AGAIN: 

:n 


THERE  ARE  3  COMPONENTS  INCLUDED  IN  RUN 

CYLINDER  TEST  COMP  B  STEEL  4.255,  5.225  CM  DIMENSIONS 

AND  THEY  COMPRISE  OF  : 

COMPONENTS 

COMP  B 

STEEL 

AIR 


ENTER  ONE  OF  THE  FOLLOWING 

CREAT  -  TO  CREATE  THE  DATE  FILE  FOR  'SIN' 

DEL  -  TO  DELETE  A  COMPONENT  FROM  THE  ABOVE 
CONT  -  TO  CONTINUE  IE.  ADD  MORE  COMPONENTS 
ALT  -  TO  ALTER  ANY  OF  THE  ABOVE  COMPONENTS 
START  -  TO  START  FROM  SCRATCH  AGAIN 
EXIT  -  TO  EXIT  THE  PROGRAM 


creat 


THE  DATA  FILE  'S INDATA '  HAS  BEEN  CREATED 


ENTER  ONE  OF  THE  FOLLOWING 


CREAT  -  TO  CREATE  THE  DATE  FILE  FOR  'SIN* 

DEL  -  TO  DELETE  A  COMPONENT  FROM  THE  ABOVE 
CONT  -  TO  CONTINUE  IE.  ADD  MORE  COMPONENTS 
ALT  -  TO  ALTER  ANY  OF  THE  ABOVE  COMPONENTS 
START  -  TO  START  FROM  SCRATCH  AGAIN 
EXIT  -  TO  EXIT  THE  PROGRAM 

•.exit 


THE  DATA  FOR  'SIN'  HAS  BEEN  PRODUCED  AND  IS  IN  YOUR  DEFAULT  ED  LIB 
' SINDATA* 

7000  job  existing 
CIO  active  mode 
:q,ed 

ED  version  15/10/82 


SIN(T400) 

GETSET ( DFC4  522) 

PURGE ( S INRUN , SN-DFC4  522,1 D-DFCEDN ) 

EXIT(U) 

REQUEST ( TAPE2 , *PF ,SN-DFC4522 

ATTACH ( SINBIN , SINBIN , SN-DFC452  2 , 1 D-DFCEDN ) 

S INBIN. 

EXIT(U) 

REWIND ( TAP E2 ) 

CATALOG ( TAPE2 , SINRUN . SN-DFC4 522. ID-DFCED IN ) 


♦EOS 

3  100  0 

CYLINDER  TEST  COM 
•20000000000E+01 
COMP  B 

50  1  0 

.84500000000E-01 
. 1 OOOOOOOOOOE-05 
.231 00000000E+00 
0. 

- . 1 2260697623E+03 
•60606060606E+00 
•30000000000E+03 
0. 

0. 

-.35222141 110E+01 
-.1 0990541 300E-01 
•65029983420E-02 
.66679850240E-01 
.lOOOOOOOOOOE+OO 
STEEL 

40  0  0 

.25000000000E-01 
.1 0000000000E-05 
.45800000000E+00 

0. 

-.1 466784021 2E+-4 
.1 26310471 14E+00 
.30000000000E+03 
.50000000000E-01 
.10000000000E+01 
AIR 

10  1  0 
.10000000000E+01 
.1 0000000000E-05 
. 3  3200000000E-01 
0. 

0. 

. 7734550235 9E+03 
.30000000000E+03 
0. 

.10000000000E+01 

0. 

0. 

0. 

0. 


0  0  0  0 
1  B  STEEL  4.225, 

.1 0000000000E+03 


.1 6900000000E-01 
. 3000000000  0E+0  3 
.1 8300000000E+01 
-.86482267660E+01 
-.341 39045850E+02 
. 5000000000 OE-04 
.1 OOOOOOOOOOE-05 
0. 

. 1 0000000000E+ 1 6 
-• 2495947472 0E+01 
-.1 543767201 OE+01 
.21 4990741 20E-03 
.1 8767838090E-01 


•50000000000E-02 
. 30000000000E+03 
.1 51 00000000E+01 
- .38238258745E+04 
-.166391 61 598E+03 
.1 1700000000E-04 
.1 OOOOOOOOOOE-05 
.1 0000000000E-05 
0. 

0 

.20000000000E+00 
. 3000000000 0E+0 3 
.1  0461 000000E+01 
.5703000000 0E+01 
0. 

0. 

.1 0000000000E-05 

0. 

0. 

0. 

0. 

0. 

0. 


10  0  0 
.225  CM  DIMENSIONS 
.40000000000E+03 


. 2  OOOOOOOOOOE+OI 

0. 

.1 OOOOOOOOOOE-OI 
-. 76497948971 E+02 
. 2  6600000000E+0 1 

0. 

0. 

0. 

.4  5000000000E+05 
.2541 1 585450E+00 
.51 978884760E+00 
.75567098200E+01 
-.1 330260321 0E-02 


.20000000000E+01 

0. 

.1 OOOOOOOOOOE-OI 
-.70321 195402E+04 
.20000000000S+01 
.70000000000E-01 
.50000000000E-02 
0. 

0. 


.20000000000E+01 

0. 

.1 OOOOOOOOOOE-OI 

0. 

0. 

0. 

0. 

0. 

0. 

0. 

0. 

0. 

0. 


0 

.20000000000E+07 


.1 6500000000E+01 

0. 

0. 

- .1 4330798590E+03 
.25300000000E+00 

0. 

0. 

0. 

.2558961 9270B-01 
.7991 8603550E-01 
- .45930275520E+00 
.50000000000E+00 


.791 70000000E+01 

0. 

0. 

-  .48367021 389E+04 
.1 0780000000E+00 
.1 500000000 0E+00 
.97800000000E+00 


0. 


.1  2929000000E-02 

0. 

0. 

0. 

0. 

0. 

0. 

0. 

0. 

0. 

0. 

0. 


10000000000E+00 


MRL-R-877 


DISTRIBUTION  LIST 


MATERIALS  RESEARCH  LABORATORIES 
DIRECTOR 

Superintendent,  Metallurgy  Division 
Dr  J.  Eadie 

Library  (2  copies) 

Mr  E.  Northeast 


DEPARTMENT  OF  DEFENCE 

Chief  Defence  Scientist/Deputy  Chief  Defence  Scientist/  (1  copy) 

Controller,  Projects  and  Analytical  Studies/ 

Superintendent,  Science  and  Technology  Programme 
Army  Scientific  Adviser 
Air  Force  Scientific  Adviser 
Navy  Scientific  Adviser 

Officer-in-Charge,  Document  Exchange  Centre  (17  copies) 

Technical  Reports  Centre,  Defence  Central  Library 
Central  Office,  Directorate  of  Quality  Assurance  -  Air  Force 
Deputy  Director  Scientific,  and  Technical  Intelligence,  Joint 
Intelligence  Organisation, 

Librarian,  Bridges  Library 

Librarian,  Engineering  Development  Establishment 

Defence  Science  Representative,  (Summary  Sheets  Only) 

Australia  High  Commission,  London. 

Counsellor  Defence  Science,  Washington  D.C.  (Summary  Sheets  Only) 

Librarian,  (Through  Officer-in-Charge),  Materials  Testing 
Laboratories,  ALEXANDRIA,  N.S.W. 

Senior  Librarian,  Aeronautical  Research  Laboratories 
Senior  Librarian,  Defence  Research  Centre  Salisbury,  S.A. 


DEPARTMENT  OF  DEFENCE  SUPPORT 
Deputy  Secretary,  DDS 

Head  of  Staff,  British  Defence  Research  &  Supply  Staff  (Aust. ) 


OTHER  FEDERAL  AND  STATE  DEPARTMENTS  AND  INSTRUMENTALITIES 


(MRL-R-877) 


DISTRIBUTION  LIST 
(Continued) 


MISCELLANEOUS  -  AUSTRALIA 


Librarian,  State  Library  of  NSW,  Sydney  NSW 

University  of  Tasmania,  Morris  Miller  Library,  Hobart,  TAS. 


MISCELLANEOUS  -  OVERSEAS 


Library  -  Exchange  Desk,  National  Bureau  of  Standards,  U.S.A. 
UK/USA/CAN/NZ  ABCA  Armies  Standardisation  Representative  (4  copies) 
Hie  Director,  Defence  Scientific  Information  &  Documentation 
Centre,  India 

Military,  Naval  and  Air  Adviser,  High  Commission  of  India, 

Canberra 

Director,  Defence  Research  Centre,  Kuala  Lumpur,  Malaysia 
Exchange  Section,  British  Library  ,  U.K. 

Periodicals  Recording  Section,  Science  Reference  Library, 

British  Library,  U.K. 

Library,  Chemical  Abstracts  Service 

INSPEC:  Acquisition  Section,  Institute  of  Electrical 
Engineers,  U.K. 

Engineering  Societies  Library,  U.S.A. 

Aeromedical  Library,  Brooks  Air  Force  Base,  Texas,  U.S.A. 

Ann  Germany  Documents  Librarian,  Hie  Centre  for  Research  Libraries , 
Chicago  II 1 . 

Defense  Attache,  Australian  Embassy,  Bangkok,  Thailand  Att.  D.  Pender 


